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We use microprobe Angle- Resolved Photoemission Spectroscopy (/iARPES) to study the Fermi 
surface and band dispersion of the Cu02 planes in the high temperature superconductor, 
YBa2Cu408- We find a strong in-plane asymmetry of the electronic structure between directions 
along a and b axes. The saddle point of the antibonding band lies at a significantly higher energy in 
the a-direction (7r,0) than the 6-direction (0,7r), whereas the bonding band displays the opposite be- 
havior. We demonstrate that the abnormal band shape is due to a strong asymmetry of the bilayer 
band splitting, likely caused by a non-trivial hybridization between the planes and chains. This 
asymmetry has an important implication for interpreting key properties of the Y-Ba-Cu-0 (YBCO) 
family, especially the superconducting gap, transport and results of inelastic neutron scattering. 

PACS numbers: 74.25. Jb, 74.72. Hs, 79.60.Bm 



It is commonly accepted that the conduction elec- 
trons/holes in the Cu02 planes play an essential role 
in determining the electronic properties of the high-T c 
superconductors. It is however not totally understood 
how the charge carriers arrange themselves and inter- 
act with each other. Strong a-b axis asymmetry of both 
the normal- and superconducting-state electronic proper- 
ties have been observed in the Y-Ba-Cu-0 (YBCO) fam- 
ily especially in YBa2Cu307_5 (Y123) and YBa2Cu408 
(Y124), which possess double Cu02 planes sandwiched 
between single and double CuO chains, respectively. 
The London penetration depth (A) measured by muon- 
spin rotation shows a strong asymmetry (A a /A& « 1.2 
in Y123±£ and A a /A 6 « 1.5 in Y124 3 ). Here the a- 
and 6-axis directions are perpendicular to and along 
the chains, respectively. Angle Resolved Photoemission 
Spectroscopy 4 observes a 50% larger gap along the chain 
direction in Y123. Angle- Resolved Electron tunneling in 
Y123/Au/Nb junctions also indicates a superconducting 
gap (A) asymmetry of A^/A a w 1.5&. Most recent re- 
ports using a similar technique, however, show only 20% 
difference 6 . Many theoretical models have been proposed 
in order to explain these results. These include the for- 
mation of striped phases 7 , orthorhombicity of the crys- 
tal structure 8 , proximity coupling between chains and 
planes 9 , broken time-reversal symmetry 10 , and an ad- 
mixture of d x 2_ y 2+s pair states 11 . The spin dynamics in 
Y123, observed in inelastic neutron scattering (INS) ex- 
periments, also shows a significant a-b axis asymmetry in 
the vicinity of the wavevector (l/2,l/2)i2ii4ii5*i£. Mook 
et. al. 13 report almost one-dimensional patterns and sug- 
gest that this is evidence for existing stripe phases 7 -. In 
contrast, recent INS results on Y123 show that the spin 



fluctuations are two-dimensional, although they have a 
strong a-b asymmetry 16 . In the latter case, a Fermi- 
liquid-based scenario might be more relevant for explain- 
ing the asymmetric spin dynamics ^ 17 ! 18 . 

Part of the controversy arises due to a lack of de- 
tailed band structure measurements in YBCO, since most 
ARPES experiments performed on this material have fo- 
cused on YBa2Cu307_5 (Y123) variety, which has frac- 
tional oxygen stoichiometry in the chains. These samples 
do not cleave well and suffer complications due to the 
ordering of oxygen in the chains. The resulting unsta- 
ble sample surfaces makes ARPES spectra difficult to 
interpret 4 . Although band calculations predict a sig- 
nificant a — b axis asymmetry of the band structure in 
Cu02 planes due to a plane-chain coupling, that is often 
ignored when interpreting the electronic properties be- 
cause of a lack of the experimental evidence supporting 
these predictions. Until now, there have been no obser- 
vation of a — b axis asymmetry in the Fermi surface of 
the Cu02 planes by angle-resolved photoemission spec- 
troscopy (ARPES) in YBCO 4 ^^ 2 ^ 2 ^ 3 ^. We chose 
YBa2Cu408 (Y124) for studying the band structure by 
ARPES, because Y124 has a fixed oxygen stoichiome- 
try that results in a much more stable surface than that 
of Y123 and avoids problems due to oxygen ordering 
in chains. We have previously demonstrated that the 
plane and chain bands can be distinguished by using 
both a small UV beam (50-100/im) and a tunable in- 
cident photon energy 24 . This lead to observation of sig- 
nificant differences in the momentum dependence of the 
bilayers splitting between Bi2212 and YBCO. Since it 
was sufficient for this purpose to acquire data only along 
nodal and a- axis directions, we could not make state- 
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FIG. 1: (Color online) ARPES intensity of a plane domain 
and corresponding EDCs along Y-S (a,b) and X-S (d,e) indi- 
cated in (g). (c,f) Symmetrized EDCs of (b) and (e) close to 
the Fermi level, (g) Fermi surfaces of antibonding-, bonding-, 
and chan-band (AB, BB, and CH, respectively) obtained by 
band calculation 27 , (h) EDCs at X and Y. (i) MDCs at the 
Fermi level along Y-S and X-S. (j) Band dispersions along Y-S 
and X-S determined from EDC peak positions of (b) and (e). 
Dotted blue lines and dashed red lines represent tight binding 
fits^ 3 - along Y-S and X-S, respectively. 



merits about the a — b axis anisotropy. In this letter, we 
compare the band structure along a and b axis with much 
improved data quality due to use of high flux beamline 
and longer acquisition times. We found a small but a sig- 
nificant a— b axis asymmetry in the band structure, which 
is mainly due to the variation of the bilayer band split- 
ting being about 40% larger along the a-direction. Our 
results suggest that the coupling between the planes and 
chains and the consequent asymmetry of the band struc- 
ture represents a crucial ingredient in our understanding 
of the origin of the asymmetries in the electronic and 
magnetic properties of YBCO systems. 

Untwinned high quality single crystals of underdoped 
YBa2Cu408 (T c ~ 80K) with a sharp superconduct- 
ing transition (AT C ~ 3K) were grown by the self-flux 
method under high oxygen pressure^. The ARPES 
experiments were carried out using a Scienta SES2002 
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FIG. 2: (Color online) ARPES intensity at momentum cuts 
parallel to Y-S (a-e) and X-S (f-j). (k-o) Band dispersion 
determined from the EDC peak positions and tight binding 
fits 33 , (p) ARPES intensity integrated within ±10meV cen- 
tered at the Fermi level as a function of k x and k y . (q) Fermi 
surface plots of antibonding-band (AB, green filled circle) and 
bonding-band (BB, purple open circle) determined from the 
peak position of the MDCs at the Fermi level. 



hemispherical analyzer mounted on the SIS beamline of 
the Swiss Light Source (SLS). As a reference for the Fermi 
energy, we used the spectral edge position of evaporated 
Au in electrical contact with the sample. The energy and 
angular resolutions were 20 meV and 0.1°, respectively. 
All spectra were measured at 25K using a 33eV photon 
energy. In order to observe both a- and 6-axis directions 
with an identical ARPES setting (photon polarization 
and geometry of analyzer and sample), we rotated the 
samples by 90° following the same micro-domain on the 
sample surface. We measured several samples and always 
obtained similar results. 

Figures 1 (a,b) and (d,e) show ARPES data (inten- 
sity maps and corresponding energy distribution curves 



3 



(EDCs)) measured along Y-S and X-S, respectively. 
Those momentum cuts are illustrated in panel (g) along 
with the Fermi surface of the planes and chains obtained 
by band calculation^. Around Y (chain direction), 
both dispersions of the antibonding (AB) and bonding 
(BB) bands, attributed to the Cu02 bilayer splitting, 
are clearly observed at higher and lower energies, re- 
spectively. Note that, due to matrix element effects, the 
ARPES intensity at positive k x is stronger than at neg- 
ative k x in the lower energy band. This k x dependence 
becomes opposite for the higher energy band. Such anti- 
correlation is typical for bonding and antibonding bands 
due to the orthogonality of their wavefunctions^. In con- 
trast to the Y direction, the AB signal is very weak close 
to X (the direction perpendicular to the chains), whereas 
the BB is clearly seen. The difference between the two 
antinodal cuts is more clearly illustrated by plotting the 
EDCs taken at X and Y in Fig.l (h), and the momentum 
distribution curves (MDCs) at the Fermi level in Fig.l 
(i). We find that the BB dispersion around X displays a 
deeper bottom (by ~40meV) and a wider Fermi closing 
than that around Y. Although a large peak correspond- 
ing to an energy state of AB is seen below the Fermi level 
(around — 20meV) in the EDC at Y (panel (h)), the peak 
edge of a weak AB signal appears to be pinned to the Ep 
at X. We used the symmetrization method 29 as well as 
the MDC peak position to determine the Fermi cross- 
ing points; EDCs are reflected about the Fermi level and 
added to the unreflected ones. This technique removes 
the Fermi function and enables us to immediately identify 
the Fermi crossing. The symmetrized EDCs of Fig. 1 (b) 
and (e) are plotted in the Fig. 1 (c) and (f), respectively. 
When the dispersion crosses Ep, two peaks in the sym- 
metrized EDC due to presence of low-lying energy states 
merge into a single one at Fermi momentum. This is 
clearly seen both for AB and BB with moving away from 
Y along Y-S (panel (c)). The superconducting gap is not 
observed in these spectra because the YBCO sample sur- 
face is known to be overdoped after cleavin g 24 ! 31 . Along 
X-S (panel (f)), the BB clearly shows a Fermi crossing. 
On the other hand, the symmetrized EDC very close to 
X has only one peak, which indicates that there is no 
Fermi crossing in AB along this cut. The shape of sym- 
metrized EDC is sensitive to the value of experimental 
energy resolution 30 . However we can confirm this from 
a behavior of the spectral intensities along T-X shown 
in Fig. 3(d); the intensity sharply decreases toward X, 
which indicates that the AB goes beyond Ep near X and 
only a small spectral tail is seen on the occupied side. 
In Fig. 1 (j), we summarize the AB and BB dispersions 
of the occupied states determined from the peak posi- 
tion of the EDCs in panels (b) and (e) along with tight 
binding fitting curves 3 - 3 -. Although our data represents 
the overdoped state, the band anisotropy reported here 
is expected to be valid for other doping levels because it 
is well known that doping of cuprates amounts to a rigid 
band shiftSa. 

Figure 2 (a-e) and (f-j) show ARPES intensity maps 
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FIG. 3: (Color online) EDCs and the symmetrized ones mea- 
sured at momenta along T-Y (a,b) and T-X (c,d) indicated in 
(h). (e) Dispersion of antibonding- and bonding-band along 
r-X (AB X and BB X ) and T-Y (AB Y and BB Y ) determined 
from peak positions of EDCs and symmetrized EDCs. Solid 
lines represent tight binding fits 33 , (f) The midpoints of AB 
and BB along T-X (Midx) and T-Y (Mid Y ). (g) Energy of 
bilayer band splitting along T-X (2t±x) and T-Y (2£xy). (h) 
Brillouin zone. 



around Y and X measured along various momentum cuts 
parallel to Y-S and X-S, respectively. Around X (panel 
(f-j)), the ARPES intensity near Ep rapidly increases 
while moving away from the antinode due to the appear- 
ance of AB on the occupied side. This is in contrast to 
the strong intensity of the AB observed for all momen- 
tum cuts around Y (panels (a-e)). We determine both 
AB and BB dispersions from the EDC peak positions, 
and plot the results near Y and X in Fig. 2(k-o); each 
panel is for two momentum cuts at equal distance from 
X and Y. The difference in the dispersions between the 
two regions reaches a maximum close to the antinodes, 
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then decreases toward the nodes. The strong a-b axis 
asymmetry of the Fermi surfaces is further visualized in 
Fig. 2(p) by plotting the ARPES intensity integrated 
within ±10meV of Ep as a function of k x and k y . In 
Fig. 2(q), we plot the Fermi crossing points extracted 
from the MDC peak positions at Ep. We find an inter- 
esting topology of the Fermi surface which appears to be 
hole-like (centered around S) close to Y and electron-like 
(centered around T) close to X. We should note that the 
previously reported data 24 was consistent with hole-like 
Fermi surface close to X. This difference is most likely 
due to small variation of doping for different batches of 
samples. In the previous work, the data was measured 
only in X quadrant and reflected about symmetry planes 
to simply show the 2D character of Fermi surface that is 
different from ID chain. In this work, we measured both 
X and Y directions on the same cleave to compare these 
in the same condition. We could get the same result for 
several samples cut out from a sample batch with high 
quality. 

In order to understand what causes the significant a — b 
asymmetry, we carefully compare the band dispersions 
along T-Y and T-X. Figure 3 (a-b) and (c-d) show the 
EDCs and the corresponding symmetrized EDCs mea- 
sured along T-Y and T-X, respectively (see panel (h)). 
We determined the energy states of BB (AB) along the 
two momentum cuts from the peak positions of the EDCs 
(symmetrized EDCs), and plot these in panel (e). We 
find that the BB along T-X has a convex downward dis- 
persion centered at X. This is in contrast with all other 
dispersions (AB along T-X and both AB and BB along 
T-Y), which show a convex upward shape. We stress 
that the former band shape is specific to Y124 because 
the latter has been observed in many cuprates^ 2 - such 
as Bi 2 Sr 2 CaCu 2 08 + 5 (B12212), Bi 2 Sr 2 Cu0 6+( 5 (Bi2201), 
and La 2 _ :E Sr :c Cu04 (LSCO). Figure 3(g) shows the split- 
ting energy (2£j_, where t± is the interlayer hopping inte- 
gral within the bi-layers) estimated along T-X and T-Y. It 
is clear that 2t± along T-X has a bump close to X and the 
difference from that along T-X is about 40% (~ 60meV). 
We therefore conclude that the characteristic band shape 



close to X (the electron-like Fermi surface in AB and a 
convex downward dispersion in BB) is a result of an en- 
hanced bilayer coupling around X, which pushes the AB 
up above Ep and pushes down the BB more than that 
around Y. In Fig. 3(f), we plot the mean-energies of AB 
and BB along T-X and T-Y. We find that the dispersion 
is almost identical along the both cuts, which indicates 
that the in-plane orthorhombicity of the Y124 crystal 
structure (b/a ~ 1.008) 25 causes only a negligible a-b 
axis asymmetry of the in — plane hopping integral (t//). 
This small crystal orthorhombicity, therefore, is unlikely 
to cause a significant a-b asymmetry of the interlayer 
hopping (£j_). The enhanced bilayer coupling around X, 
where the Fermi surface of chains is adjacent to that of 
planes, is most likely due to a non trivial (momentum de- 
pendent) coupling between the chains and plane a 26 i 27 . It 
is worth noting that the current value of 2t± (~160meV) 
is about 4 times smaller than that estimated from band 
calculations (~600 meV)21, which indicates that strong 
electron correlations suppress the interlayer hopping. 

In conclusion, we find a strong in-plane a — b asym- 
metry of the electronic structure in YBa 2 Cu408 (Y124). 
The Fermi surface in the antibonding band has a hole-like 
shape (centered around S) close to Y and an electron-like 
shape (centered around T) close to X. This asymmetry 
is caused by an enhanced bilayer band splitting close to 
X, where the Fermi surface of chains is adjacent to that 
of planes, likely due to a nontrivial coupling between the 
chains and planes. These results are crucial to under- 
stand the abnormal asymmetries observed in the super- 
conducting gap and the significantly anisotropic spin dy- 
namics, strongly supporting a number of recent theoret- 
ical work G 8 i 17 i 18 . 
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We use seven tight-binding fitting functions with a — b 
axis anisotropy: no = 1, n\ = [(1 + Si) cosk x + cos/c y ]/2, 
rt2 = cos k x cos k y , 713 = [(1 + Ss)cos2k x + cos2A; y ]/2, 
714 = [(1 + S4) cos 2k x cos k y + cos k x cos 2k y ]/2, 
n$ = cos 2k x cos 2k y , uq = [(1 + 5q) cos3k x + cos3/c y ]/2. 
The fitting parameter (to, ti, ti, £3, £4, £5, ^6, Si, S3, £4, Se) 
for obtained tight-binding band [e(k x ,k y ) — ^tira] 
is (0.242,-0.582,0.189,-0.055,-0.023,-0.007,-0.019,0.029,- 
0.073,-0.300,0.334) and (0.100,-0.665,0.166,-0.102,-0.108,- 
0.009,-0.055,-0.002,-0.072,0.350,-0.056) for antibonding- 
band and bonding-band, respectively. 



